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1. Introduction and Summary 

Recently the new very interesting formulation of the non-linear massive gravity |2| was 
introduced with significant improvement reached in ||3|, |^. This theory was further extended 
in where the theory was formulated with general reference metric. The most crucial 
fact that is related to given theory is the proof of the absence of the ghosts that are 
generally expected in any theory that breaks the diffeomorphism invariance. As is well know 
the physical content is determined in the Hamiltonian formulation when all constraints 
are identified together with their nature. This analysis was performed in several papers 
||6|, 1^, ^, ^, 10, 11, 12] with the most important results derived in ]T^ , p^ ] with the outcome 



that this non-linear massive theory possesses one additional constraint and the resulting 
constraint structure is sufficient for the elimination of the ghost degree of freedom. 



Very interesting extension of given theory was suggested in |15] when the kinetic term 



for the general reference metric was introduced and hence g^j^ and /^^ come in the sym- 



metric way in the action. Then it was argued in |13] that the resulting theory is the ghost 
free formulation of the bimetric theory of gravity ^ . 

The goal of this paper is to perform the Hamiltonian analysis of the bimetric theory 
of gravity in the form introduced in ||l5|] . For simplicity we call this theory as the new 
bimetric theory of gravity (NBTG). We would like to explicitly determine the structure 
of the constraints and eventually to prove the absence of the ghosts. Remarkably we find 
very subtle issue related to NBTG which forces us to doubt whether the ghost could be 
eliminated in given theory or not. More explicitly, the non- linear massive gravity with 
general reference metric has the potential that depends on the matrix H^,/ = g^^fpcr where 
fpa is fixed background metric. There is now no doubt that such theory is ghost free. 
On the other hand the situation changes in case of the bimetric theory of gravity when we 
promote f^i, as an additional dynamical field with the kinetic term given by Einstein-Hilbert 
action. Then, since the interaction term between two metrics has the square structure 



^For further analysis of given theory, see 0, |l|, ||, ||, §||4[ H H, ||, §|. 
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form we follow |13, 14, ^ and perform the redefinition of one shift function that makes 
the theory linear in A'^ and M which are the lapse functions in and /^j^ respectively. 
It is important that the action has the same structure as the non-linear massive gravity 
action with general reference metric f^y with additional kinetic term for /^j^. However the 
fact that ffj_i, is dynamical has crucial impact on the Hamiltonian structure of given theory. 
Explicitly, the Hamiltonian is given as the linear combination of the constraints as opposite 
to the case of the non-linear massive gravity where the Hamiltonian does not vanish on 
the constraint surface. Now the crucial point is that the components of the metric /^jy 
that appear as the fixed parameters in the non-linear massive gravity case should be now 
considered as Lagrange multiplicators whose values are determined by the requirements of 
the preservation of all constraints during the time evolution of the system. However then we 
find that the requirement of the preservation of the constraint Cq whose explicit definition 
will be given below leads to the differential equation for the Lagrangian multiplicator M 
that is related to the P. In other words the value of the Lagrange multiplicator M is 
determined by the requirement of the preservation of Cq during the time evolution of the 
system. In the same way we fix the value of the Lagrange multiplicator A^. Say differently, 
Co and T> could be interpreted as the second class constraints. 

At this place we should compare our result with the known proof of the absence of 
the ghosts in the bimetric theory. It was shown in the very nice paper that the 
requirement of the preservation of the constraint Cq implies an additional constraint C(2) 
(in their notation) given in the e.q. (3.32) in this paper. We see that this constraint 
contains the covariant derivative of M which, as we argued above, is fixed in case of the 
non-linear massive gravity theory so that it is really natural to interpret C(2) additional 
constraint. Then this constraint together with Cq are responsible for the elimination of the 
ghost mode which is crucial for the consistency of non-linear massive gravity. However in 
case of the bimetric gravity M should be considered as the Lagrange multiplicator whose 
value is fixed by the consistency of given theory. 

The fact that Cq and D should be interpreted as the second class constraints has im- 
portant consequence for the dynamics of the theory. More precisely, since the Hamiltonian 
is given as the linear combination of the constraints implies that the resulting Hamiltonian 
vanishes strongly. This is rather puzzling result and we believe that this is a consequence of 



the redefinition of the shift function [03, 14, 15] which is certainly useful for the non-linear 



massive gravity where the diffeomorphism invariance is explicitly broken but we are not 
sure whether it is suitable for the bimetric theory of gravity where all fields are dynamical 
and the theory is manifestly diffeomorphism invariant under diagonal diffeomorphism. In 
fact, it is non-trivial task to identify such generator as was shown recently in [^] in case 



of particular model of bimetric theory of gravity |35|. We are currently analyzing NBTG 
following [^] and we believe that it is possible to identify four first class constraints that 
are generators of the diagonal diffeomorphism. On the other hand the analysis performed 
so far suggests that it is very difficult or even impossible to find an additional constraint 
that could eliminate the additional mode 



^This analysis will appear in forthcoming publication. 
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We should however stress that we have to be very careful with definitive conclusions. 
We wanted to show that the extension of the non-linear massive gravity to the bimetric 



theory of gravity as was performed in |15] could be more subtle that we initially thought 
and that it is not really clear whether given theory is ghost free or not. It is still possible 
that there exists the way how to find four first class constraints that are generators of 
diagonal diffeomorphism together with additional constraints that eliminate ghost mode. 
Clearly more work is needed in order to resolve this issue. 

This paper is organized as follows. In the next section (pi) we review the Hamiltonian 



formulation of the bimetric theory of gravity [15| and identify primary and the secondary 
constraints. Then in section (|3|) we calculate the algebra of constraints for the case of 
the minimal version of the bimetric gravity. Finally in Appendix ^ we briefly discuss the 
Hamiltonian formulation of the bimetric F{R) theory of gravity. 

2. Bimetric Gravity 

In this section we review the main properties of the bimetric theory of gravity in the 
formulation presented in [15|. The starting point is following action 



where 



n=0 



1 1 



(2.1) 



and where gf^,y, f^i, are four-dimensional metric components with ^'^^R^^\ ^'^^R^^^ correspond- 
ing scalar curvatures. Further, efc(A) are elementary symmetric polynomials of the eigen- 
values of A. For generic 4x4 matrix they are given by 

eo(A) = 1 , 
ei(A) = [A] , 

e2(A) = i([A]2 - [A^]) , 

e3(A) = i([Af -3[A][A2] + 2[A3]) , 

e4(A) = ^([A]4-6[Ap[A2] + 3[A2]2 + 8[A][A3]-6[A4]) , 

efc(A) = ,for A; > 4 , (2.3) 

where A^u is 4 x 4 matrix and where 

[A] = TrA = A^^ . (2.4) 
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Of the four /3„ two combinations are related to the mass and the cosmological constant 
while the remaining two combinations are free parameters. If we consider the case when the 
cosmological constant is zero and the parameter m is mass, the four ^„ are parameterized 
in terms of the as and 04 of Q 

/3„ = (-l)"Q(4-n)(3-n)-(4-n)Q3 + a4) • (2.5) 

The minimal action corresponds to /32 = /^s = that implies 03 = 04 = 1 and consequently 
/3o = 3 ,/3i = -1. 

Our goal is to find the Hamiltonian formulation of given theory and determine cor- 
responding primary and the secondary constraints. As the first step we introduce 3 + 1 
decomposition of both g^j^i, and f^i, ^] 

500 = -N"^ + Nig'^Nj , goi = Ni , gij = gij , 



-00 _ _J_ ^Oi ^ ^ -ij ^ ij _ N'N^ 



(2.6) 



and 



/oo = -M^ + Lif^Lj , foi = Li , fij = ft 



rOO _ fOi _ fij fij Ti T fji 



(2.7) 



and where we defined g"^^ and f^^ as the inverse to gij and fij respectively 

= , f,uf^^ = 5/ . (2.8) 

Following ^ 1^, ^, ^ we perform following redefinition of the shift function 

A^^ = Mh' + V + ND'jh^ (2.9) 

so that the resulting action is linear in M and A^. Note that the matrix D^j obeys the 
equation |, |, |, |ll 



V^&j = yj{g'^ - &„,n"^D\h-)fkj (2.10) 
and also following important property 

fikD^j = fjkD\ . 

(2.11) 

Then after some calculations we derive the bimetric gravity action in the form |^ ^ ^, |15[ 
S = M} j dtSi^M^f[K,,g'^''^kki + R^f^ + Ml j dtd'^^N^lKijG'^^^Kki + iJ^^^] + 
+ 2m'^Mlff j dtd^x^{MU + NV) , 

(2.12) 
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where 

Kij = -^{dtgij - ViNj{h,g) - VjNi{fi,g)) , 

(2.13) 

where 

Ni = Mgijh^ + gijU + NgikD^'jU^ , U = fijL^ , (2.14) 

and where Vi,i?(s) and Vi,R^f^ are the covariant derivatives and scalar curvatures calcu- 
lated using gij and fij respectively. Further, Q'^^^'- and Q^ki g^^g ^jg Witt metrics defined 
as 

with inverse 

Qijki = -{gikgji + gugjk) - -gijgu , Qijki = -{fikfji + fufjk) - -fijfki (2-16) 

that obey the relation 

QijkiQ = -{Oi dj +di dj ) , GijkiG = -{di dj +di dj ) . (2.17) 
Finally, V and U introduced in (|2.12|) have the form 



(2.18) 

where 

X = l-h'fijh^ . (2.19) 

The action ( |2.12| ) is suitable for the Hamiltonian formalism. First we find the momenta 
conjugate to A^, and gij 

7r,v~0, TTi^O, t:'^ = Ml^g'^'''Kki (2.20) 

together with the momenta conjugate to M, U and fij 

PM-0, p,^0, p'^ =M}^-g'i^'kki . (2.21) 
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Then after some calculations we find following Hamiltonian 

= J d^x(iVCo + MV + VTZi) , 



(2.22) 



where 



(2.23) 

where we also denoted 

7^S^) = -2g,kVy^ , = -2f,kVip"' . (2.24) 

From previous analysis we see that we have eight primary constraints 

TTAT « , TT j « , pM ~ , Pi^O . (2.25) 

Then the next step is to analyze the requirement that these constraints are preserved during 
the time evolution of the system 

dtT^N = {VTAT, H} = -Co « , 
dtpM = {PM, H} = -T> « , 

dtTTi = {7ri,H} = Ck MSf + AT \^0, 



where 



dtPi = {Pi,H} = -TZi « 



(2.26) 



'2 



(2.27) 



and where we used the canonical Poisson brackets 

{iV(x),vr^(y)} = 5(x-y) , {n^x), 7r,(y)} = (5}5(x - y) 

[g,^{^U^\y)] = l(5f4+5^5^)5(x-y) , 
{M(x),pM(y)} = <5(x-y) , {LXx),p,(y)} =4<5(x-y) 

{/,,(x),p'^'(y)} = ^-{6^s\ + 6l6^)6{^-y) 



(2.28) 
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and also following important relations |13] 



5^D\ 1 



c--,- — fnm i 

orv yjx Orv 



(2.29) 



In summary we have following 16 constraints 



primary : ttat « , TTj « , /o^i/ , pi , 
secondary : Cq ^ , P , Cj , T^j f« . 



(2.30) 



Now we have to check the stability of all constraints when the total Hamiltonian takes the 
form 

Ht = j d^x{NCo + MV + L'TZi + u^ttn + u'th + 

(2.31) 

where A^, M, L*, n^, n*, f v*, S* are Lagrange multiplicators related to the constraints 
(2.30). For simplicity we restrict ourselves to the case of the minimal bi-metric theory. 



3. Preservation of Constraints in Case of Minimal Bimetric Gravity 

The minimal bimetric theory is defined by the following choice of parameters 

/3o = 3,/3i = -l,/32 = 0,/33 = 0,/34 = l. (3.1) 



Now we proceed to the analysis of the preservation of all constraints given in ( |2.3C1| ). It 
is easy to see that the constraint vr^v ~ is trivial preserved. On the other hand the 
requirement of the preservation of the constraint tTj ~ takes the form 



M6f + 



dfi'- 



Ck+ / d^X^^ {TTiXji^)} =0 



(3.2) 



where 



{7ri(x),Cj(y)} 



5{y.-y) = A^^,c,5(x - y) . 



(3.3) 
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Since 

det A^^,c, = det (^^^ det (^6) + in'^n'"/™,^ = ^ det fi, / 

(3.4) 

we find that A,r,;,Cj is non-singular matrix on the whole phase space. However this fact also 
implies that the equation ( p.2| ) has trivial solution 

= . (3.5) 

In the samy we proceed with the analysis of the time evolution of the constraint Cj 

dtdi^) = {C,{x),Ht} = I dV(iV(y){C.(x),Co(y)} + 

+M(y){C,(x),P(y)}+t;^(y) {C,(x), 7r,(y)}) = . 

(3.6) 



According to ( |3.4| ) we see that ( |3.6[ ) can be solved for as functions of the canonical 
variables and N, M. Say differently, vrj and C, are the second class constraints. 

As the next step we consider the constraint TZi. It turns out that is convenient to 
extend it by the expression dih^iTj + dj{h^iTi) and consider its smeared form 

Ts{N') = J cfxN\n^^ + 7^^^) +p^di(^ + dih^TTj + dj{fi^TTi)) = J cf^N'TZi . (3.7) 

Then using the canonical Poisson brackets we find 

{Ts{N'),g^j{^)} = -afc5ii(x)iV'=(x) - a,iV'=(x)r7fc,(x) -5ifc(x)a,iV^(x) , 

{T5(iV*),Vr^^'(x)} = -dk{N\^)7T'^^))+dkN\^)7T'^{^)+7T''{^)dkN^{^) , 

{T5(iV*),p^^'(x)} = -dk{N'^{^)p^^{^)) + d,N'{^)p'^{^)+p^H^)dkN^(.^) , 
{Ts(7V^),n*(x)} = -N''{x)dkn'{x) + djN\x)h^x) , 
{Ts(iV*),^i(x)} = -a,.(iV'=^i)(x)-aiAr^Xx)vrfc(x) . 

(3.8) 

Then we easily find the familiar result 

{Ts{N'),Ts{M^)} = TsiiN^djM' - M^djN')) . (3.9) 

To proceed further we need to know the Poisson bracket between Tsi^'') and D^j which 
can be determined when we know the explicit form of D'^j ^ ^, ^, |l5| 



(3.10) 
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where 

Q\ = x5'^ + n^nVfc, , {Q-^Y t = h^'k " n'n'^fmk) • (3.11) 



Using the exphcit form of ■ given in (3.10) we see that the Poisson bracket between 



Ts(A^*) and • is determined by the Poisson brackets between T5(A^*) and gij,fij and 



The Poisson brackets between T5'(A^*) and gij and fij were given in (|3.8D and the 
Poisson bracket between Ts(A^*) and Q^j can be easily determined using ( ^ and (|3lll ) 



{Tsm,Q'j} = -dkQ'^N^ + dkN'h^n^fm.-n'n'^fmkdjN^ = 

= -dkQ'jN'' + dkN'Q'^j - Q\djN'' . (3.12) 

Then with the help of this result we find 

{t5(7V^), Z)^} = -dk&jN^ + dkN'D) - D\djN^ (3.13) 

and finally collecting all these results we obtain 

{Ts(iV*),Co} = -diCoN' - d.N'C^ , 
{Ts{N'),V} = -diVN' - d^N'V , 
{TsiN'),Ci} = -djN^d - N^djC^ - diWCj . 

(3.14) 

Then it is easy to see that TsiN"^) is preserved during the time evolution of the system 
and that it corresponds to the generator of the spatial diffeomorphism. In other words TZi 
are first class constraints. 

Now we come to the calculation of the Poisson brackets between the constraints Co 
and P. It turns out that it is useful to introduce the smeared form of these constraints 

C{N) = j d^xiV(x)Co(x) , D(M) = j d^xM{x)V{x) . (3.15) 

We begin with the Poisson bracket between Co(x) and Co(y). Since Cq does not depend on 
/?*-^ we immediately find that the Poisson bracket between Co(x) and Co(y) has the same 
form as in |13] which means that it vanishes on the constraint surface 



{Co(x),Co(y)} «0 . (3.16) 

In case of P we find 
{D(M),D(iV)} = 

= j (fx{^^MN - diNM)[f ^nj + {h'h^ - f^)Cj] . 

(3.17) 
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We see that the right side vanishes on the constraint surface. Finally we come to the 
calculation of the Poisson bracket between C(M) and D(A^) 

/[' ( Am'^M'^ 



where 



+ J d^:!i[N&^h'^h'n[^^djM - Mh^^jND\^fl"'^zl^^]- 
- ^m^Mlfj j dMNV'''Vi{Mh%,k - Vp{Nb\n')gkmU'^PM\ 



(3.18) 



^^kl ' X fm.n ' V J 



§gkl ' ^fmn ' ^""dgkl ' 

(3.19) 

Let us analyze the Poisson bracket calculated above in more details. First of all we see 
that the first expression vanishes on the constraint surface Cq ~ which is desired result. 
On the other hand in order to analyze the time evolution of the local constraint Co it is 
useful to express the local form of the Poisson bracket from ( 3.18| ) that can be schematically 
written as 

{C(iV),D(M)} = 

d^z(7V(z)M(z)F(z) +a^.iV(z)V*(z)M(z) + A^(z)9^.M(z)W*(z)) , 

(3.20) 

where the explicit form of F, V*, W* follow from ( p.l8| ). Let us now write 

N{z) = J d'^xN{x)6{x - z) ,M(z) = J d^yM{y)6{y - z) (3.21) 

and insert it to the right side of the Poisson bracket ( p. 20 ). Then after some calculation 
we find that it is equal to 

J dWyiV(x)M(y)[5(x - y)F(x) + ^^(x - y)V(y) + ^^(x - y)W*(x)] . 

(3.22) 

On the other hand we have 

{C(iV),D(M)} = J d3xdViV(x)M(y){C(x),2?(y)} . (3.23) 
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Since ( 3.22| ) and ( |3.23| ) have to match for any A^(x),Af(y) we obtain 



{C(x),P(y)} = 5(x - y)F(x) + ^^(x - y)V(y) + ^<^(x - y)W^(x) . (3.24) 

Now using this expression we can easily determine the requirement of the preservation of 
the constraint Cq during the time evolution of the system 

aA(x) = {Co{^),Ht} « I d-VM(y){C(x),P(y)} = 

= M(x)[F(x) - 5,VXx)] + ^^[VF^x) - y^x)] = 

(3.25) 

This is the most crucial point of our calculation that deserves careful explanation. Let us 



imagine that we have V* = W*. Then ( p.25D has solution either M(x) = or -F(x) — 
9iV*(x) = 0. In fact, the first case occurs when the expression -F(x) — 5jV*(x) is non-zero 
on the whole phase space, as for example in case when this expression is constant. On the 
other hand when -F(x) — 5jV*(x) depends on the phase space variables it is more natural 
to impose the condition Cq^^ = -F(x) — c?iV*(x) = as an additional constraint. This 
would be the desired result since now we would have two second class constraints Cq,Cq^ 
that would be sufficient for elimination of the ghost mode. Unfortunately as we can see 
from ( 3.18| ) V* 7^ W* and the situation is completely different since the equation ( 3.25| ) 



cannot leave M undetermined. Rather we should interpret ( p.25| ) as equation that can be 
solved for M as function of the phase space variables. In fact, in the same way we can 
analyze the requirement of the preservation of the constraint T) that again leads to the 
differential equation that can be solved for A'^. In other words we mean that it is now 
natural to interpret Cq together with T> as the second class constraints. Certainly this is 
very strange result. In particular, now we find that the total Hamiltonian strongly vanishes 
up the diffeomorphism constraint. Of course, we know that this cannot be right since the 
theory possesses the overall diffeomorphism invariance and hence there should be four the 
first call constraints that are generators of this diffeomorphism. The way how to find 
such generators for bimetric theory of gravity was suggested in [|l^] at least for particular 
bimetric gravity model. The extension of this work to the case of the non-linear bimetric 
gravity is currently under consideration. Then the result derived in this section suggests 
that the redefinition of the shift function which is very useful in the case of the non-linear 
massive gravity may not be the right way in the case of the bimetric theory of gravity. 

Despite of the fact that the total Hamiltonian vanishes it is instructive to count 
the number of the physical degrees of freedom. Recall that phase space variables are 
A^, tttv, ra*, vTj, M, pm, L*, Pi, /jj, p*-', ^fjj, vr*-' so that the total number of the phase space de- 
grees of freedom is Np^s.d.f. = 40. On the other hand we have iV/.c. = 8 first class constraints 
vr^v ~ , pm ~ , Pi w , T^j 0. Finally we have Ns.c. = 8 second class constraints 
Co ~0 ssO ,7rj ssO ,Cj ssO. Then the number of the physical degrees of freedom is 



equal to |32] 



Nf.d.f. = Np,s.d.f. - 2Nf,^, - Ns.c. = 16 . (3.26) 
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At linearized level we can identify four degrees of freedom corresponding to the massless 
graviton, ten degrees of freedom corresponding to the massive graviton and two additional 
degrees of freedom corresponding to the ghost mode. It is important to stress that the 
same result can be found when we identify four first class constraints corresponding to the 
diagonal diffeomorphism and also additional eight second class constraints as in case of 
the bimetric gravity model analyzed in [19|. Of course, the square root structure of the 



potential has remarkable property in case of the non-linear massive gravity and maybe it 
could be useful in case of the bimetric gravity as well. We only say that the step from the 
non-linear massive gravity to the bimetric gravity is not straightforward as it seems to be. 

A. Hamiltonian Analysis of F{R) Bimetric Gravity 

In this appendix we briefly perform the Hamiltonian formulation of F{R) bimetric theory 



of gravity which was introduced by S.Odintsov and Nojiri in [22|. 



The starting point is the action for the non-linear bimetric gravity theory 



n=0 

(A.l) 

Then in order to construct the F{R) analogue of the bimetric massive gravity we add 



following expression to the action (A.l) 



5i = -Ml I d^x^g {^g^''d^ct>d,<l^ + 1/(0)) . (A.2) 



Then with the help of the Weyl transformation 



R[g] = e*{R[g'] - \g'^^V'^4>V',^ + S^'^'^V^^V^) 



(A.3) 



we find that Stot = -|- 5*1 takes the form 

SpK = M] I d^x^fR^f) + Ml I d^x^[e-'f'R[g'] - e-^'f'V{cP)] + 



2m'M^jf I d4xv^^/3„e(t-2)<^e„(yF^) 

n=0 



using 



(A.4) 

.(\/rV) = e't>l^en{4^f) ■ (A.5) 
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In what follows we will consider ( [A. 4 ) as the definition of the F(R) bimetric theory of 
gravity. Of course we should be able to solve the equation of motion for (j) at least in 
principle so that we could express (f) as the function of g' , f ■ Then inserting back to the 
action ( A.4D we obtain the action that is non-linear function of R[g'] and hence has the 
form of the F(R) theory of gravity ^. For that reason we can name ( [A .41 ) as bimetric F{R) 
theory of gravity even if its definition using the scalar field is more natural. Finally, in the 
following we omit ' over g^^, ffiv 

To proceed further we perform the redefinition of the shift A^* as in section (^) so that 
we find the F[R) bigravity action in the form 

Sfr = Mj J dtcf^lM^K.^g'^'^^Kki + v^MRf] + 

+ m2 J dt(f^^N[e-^K,jG'^^^Kki + e-'^R^3) _ e-^'^V^ + 

+ 2m^Mljf j dtcfx^N{MU + NV) . 

(A.6) 

Note that ( A. 61 ) has similar form as the action (2.12) up to presence of the additional scalar 
potential V{cj)) and powers of the factor e*^. Explicitly, we have 

U = /3ie-t*V5 + P2e-*Wfb\ + h'Ujb\h^] + 

+ P^e-\'t>Wi{b\n'h,b\n^ - b\n^ U,b\n') + ]^^\b\br - & ^b\)] + (3,^ . 

(A.7) 

Now using the action ( |A.6| ) we can find the corresponding Hamiltonian. Firstly we find 
the momenta conjugate to A^, and gij 

vr^^O, vr.^O, t:'^ = M^^e-'t'O'^^' K^i (A.8) 

together with the momenta conjugate to N, Li and fij 

PM^O, p'^ =M}^fg'^^'ku . (A.9) 

Since the action (A. 4) does not contain the time derivative of (j) we find that the momentum 
conjugate to 4> is zero 

P4,^Q . (A. 10) 

As a result we find following Hamiltonian 



H = j d^^{7T'^dtg.,j + p'Wthj - C) 



^For review, see 133, 34 
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d^x(ArCo + MV + L'TZi) , 

(A.11) 

where 



(A.12) 



Comparing with the situation in the second section we see that there is an additional 
primary constraint ptf,^ 0. Then again the requirement of the preservation of the primary 
constraints imphes the secondary constraints that have the same form as in case of pure 
bimetric theory of gravity. There is however an additional constraint Q that follows from 
the requirement of the preservation of the constraint « 

dtp^ = {p^,H} = N (^--^T,^^g^^^,^i + e-^^M^R^<^^- 



(A.13) 



In summary we have following set of 18 constraints 



primary : ttjv ~ , TTj , pM ~ , , 
secondary :Co~0,PRiO,CiRiO,7?.iRiO,^RiO. 



Note that now the constraint Ci has explicit form 

d = +2m^M^ff^x 



(A.14) 



(A.15) 



Now we should check the stability of all constraints when the total Hamiltonian takes the 
form 

Ht = J d^x{NCo + MT> + VTli + u^p^ + -u^ttat + nVj + 

+ v^pM + v'pi + u'^'g + s^Ci) . 

(A.16) 
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It is easy to see that the Hamiltonian structure of F{R) bimetric theory of gravity is 
almost the same as the structure of NBTG analyzed in previous two sections with small 
exception that there are two additional constraints ^ ,Q 0. They are the second 
class constraints that vanish strongly and can be explicitly solved with respect to and 
(p at least in principle. On the other hand they do not affect the analysis of all remaining 
constraints so that the constraint structure of given theory is the same as in case of non- 
linear bimetric gravity. For that reason we will not repeat the calculations performed in 
section (^). 
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